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y Based on an accurate computation of the first 1851 quantal energy levels of the 

truncated hyperbola billiard, we have found an anomalous long-range modulation in 
(— | the integrated level density. It is shown that the observed anomaly can be explained 

t~) by an additional term in Gutzwiller's trace formula. This term is given as a sum over 

families of closed, non-periodic orbits which are reflected in a point of the billiard 
boundary where the boundary is continuously differentiable, but its curvature radius 
changes discontinuously. 
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During the last few years a lot of work has been devoted to the semiclassical 
quantization of classically chaotic Hamiltonian systems - especially Euclidean bil- 
liards and the geodesic flow on negatively curved Riemannian surfaces. The main 
effort of these examinations lay in probing Gutzwiller's trace formula, see eq. (8), 
that was hrst derived in [1]. For a classically chaotic billiard system with a smooth 
boundary and unstable, isolated periodic orbits only, this formula gives a semiclas- 
sical relation between geometrical properties of the periodic orbits and the quantum 
mechanical energy spectrum. For certain surfaces of constant negative curvature, 
this formula coincides with the exact, rigorously proven Selberg trace formula of 
analytic number theory [2]. 

Recently billiard systems have been examined which break the conditions stated 
above for applying Gutzwiller's trace formula. One example is the well-known Buni- 
movich stadium billiard [3]. In this system there exists a family of neutral peri- 
odic orbits (so-called bouncing ball orbits), and in addition the boundary curvature 
changes discontinuously in four points, which we will denote by critical points. In 
[4] and [5] the contribution of the bouncing ball orbits to the quantum mechanical 
level density has been calculated. Also, the effect of the boundary discontinuities 
has been studied in these papers. In the stadium billiard it turns out that periodic 
orbits, which are reflected in one of the critical points, need special consideration. 
In [4] and [5] the contribution of their single traversal has been calculated for the 
desymmetrized and the full stadium billiard respectively. In this Letter we derive the 
semiclassical contribution of non-periodic "critical orbits" to the trace formula and 
show that under certain conditions their contribution even exceeds the periodic-orbit 
contribution. 

The investigated system is the truncated hyperbola billiard, see fig. I, which is a 
variation of the unbounded hyperbola billiard studied previously in [6]. Its boundary 
consists of the x-axis, the diagonal y = x, and the hyperbola y = -, truncated at the 
point A = (a, -) by a circular arc perpendicular to the x-axis. At the connection 
point A, the boundary is continuously differentiable, but its curvature radius changes 
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discontinuously from R + = — f f^ 2 + ;^) 2 (hyperbola) to R~ = (^2 + ^e) 2 (circu- 
lar arc). In our numerical calculations presented below, the value of the parameter 
a is chosen to be a = 7.1 1 75806. The motivation for investigating this billiard (with 
the given value for a) stems from a recent experiment [7] with a flat superconducting 
microwave cavity shaped like the truncated hyperbola billiard. For frequencies be- 
low approximately 20 GHz, only TM modes are excited which are described by the 
two-dimensional Helmholtz equation —AE Z = k 2 E z with Dirichlet boundary condi- 
tions. This wave equation coincides with the two-dimensional Schrodinger equation 




Figure I: Shape of the truncated hyperbola billiard for a ss 7.12 



— ^Atjj = Eip for the quantized billiard system. A detailed comparison of the 
measured resonances with the computed eigenvalues will be published elsewhere. 

The Green function for the Schrodinger equation with Dirichlet or Neumann 
boundary conditions can be approximated semiclassically ifi — > 0) by 

G{q",q',E) Pd Y, T=V^ 

cl. orbits iW2nih 

Xexp^S(q",q',E)-i^fiY (1) 

This relation can be derived from the well-known Van Vleck formula by Fourier 
transformation [1]. From this one obtains the semiclassical contribution of classi- 
cal orbits to the quantum mechanical eigenvalue spectrum by taking the trace (0 
denotes the billiard domain) 

g(E)= I cPqG(q,q,E), (2) 

Jn 

whose imaginary part is proportional to the quantum mechanical level density 

oo -i 

d(E) = J2${E- E n ) = —lmg(E). (3) 

According to eq. (1), the Green function is approximated by a sum over all classical 
trajectories running from q' to q" with energy E. The phase of these oscillating 
terms is given by the classical action S(q" } q' } E), and an index // counting twice the 
number of reflections on Dirichlet boundaries plus the number of conjugate points 
along the classical orbit. The amplitude is given by the square root of Van Vleck's 
determinant 

' d 2 S/(dq'dq") d 2 S/(dEdq") ' 
d 2 S/(dq'dE) d 2 S/dE 2 

which in two dimensions can be expressed by the upper right element of the classical 
2x2 monodromy matrix M: 



» = -**{ V'Z-,7, !S' . W 



J_^5_ = J_J_ 

v'v" dq' 2 dq'l v'v" M 12 ' V ; 

Here v' and v" denote the particle's velocity in q' and q" respectively. The definition 
of the monodromy matrix and a method of computing it in case of billiard systems 
can be found in [6]. For billiard systems with twice continuously differentiable 
boundary, Af i2 is a continuous function of energy and the points q' and q" . The 
integration over configuration space in eq. (2) can therefore be approximated by 
focusing one's attention on the neighborhood of points q = q' = q" with a stationary 
phase, i.e. 



dS(q",q',E) | dS(q",q',E) 
dq" dq' 



f - p> = 0. (6) 
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Figure 2: Example of a critical orbit with starting point P, and the critical curve 
(broken line) of its family 

Thus the stationary points are the initial and final positions of trajectories which are 
closed in phase space - the periodic orbits. In their neighborhood, the slowly varying 
amplitude \J\D\ can be regarded as constant and the action can be expanded up to 
second order in the coordinate variable q 2 which is perpendicular to the orbit itself: 

S(q,q,E) = S+^Kq 2 2 + --- (7) 

Following this program, one eventually obtains Gutzwiller's trace formula [I] 
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where J2 po 1S a sum over all periodic orbits with energy E, L is the geometrical 
length of the primitive periodic orbit (po), and a = ji — |(sgn(/c) — I) is a coordinate 
invariant Morse index attached to the periodic orbit [8]. 

The boundary of the truncated hyperbola billiard contains a so-called critical 
point at A = (a, -), where the boundary curvature changes discontinuously. As the 
monodromy matrix of a classical trajectory depends on the boundary's curvature 
in the reflection points, the closed, but non-periodic orbits being reflected in such 
a critical point have to be examined thoroughly. We will call them critical orbits 
(co), see fig. 2. They are typically forming one-parameter families. This can be 
seen by continuously varying the angle of incidence at the critical reflection point 



A. The set of initial and final positions of such a family F of closed critical orbits 
yields a smooth curve Cf which we call critical curve, see fig. 2. In their vicinity 
the stationary phase method fails due to the discontinuous amplitude factor J\D\. 

We now determine the contribution of critical orbits (strictly speaking: of orbits 
in their vicinity) to the trace g(E). This can be done analogously to the derivation of 
Gutzwiller's trace formula indicated above. A curvilinear coordinate system (si, s 2 ) 
is introduced along each critical curve, si being the arc parameter and s 2 the coor- 
dinate perpendicular to the critical curve. Along the direction of s 2 , two separate 
integrations have to be performed - for the orbits with the crucial reflection point to 
the left and to the right of the critical point A respectively. In each ^-integration, 
the matrix element M\ 2 is a continuous, slowly changing function of s 2 and can thus 
be approximated by the constant value M 12 = lim S2 ^ Afi 2 with ± denoting the 
different values for the two families of orbits, depending on the different curvature 
radii R + and R~ in the critical point. 

Now the action is expanded up to second order in s 2 

S(s, s, E) = S + 0s 2 + ^Ks 2 2 + ■■■ (9) 

(note that 6 ^ for non-periodic orbits), and the integration interval is extended to 
[0, oo). At this point, special care has to be taken in the case of non-periodic orbits. 
For sgn(#) ^ sgn(«) there is a stationary point of the approximated classical action 
at the point s 2 = — inside the integration interval. But either the genuine action 
has no stationary point at all, or its contribution is already included in Gutzwiller's 
trace formula (8) for periodic orbits. In both cases the integral J^° ds 2 has to be 
regularized by subtracting from it the contribution ff^ ds 2 of the "false" stationary 
point. 

Keeping this in mind, the ^-integrations can be performed explicitly resulting 
in the sum of two error functions. The details of these calculations can be found in 
[9]. The final result is 

yJc F Sz/^ti^M^I 

f % 7T 

xsgn(/€ e ) exp I —S — i—a e — iU ( 



xeifc (V-iU e ) . (10) 

In this trace formula, J2f 1S a sum over all one-parameter families F of non-periodic 
critical orbits. f c ds\ is taken along each critical curve and cannot be evaluated 
analytically (in contrast to the periodic-orbit theory). It is therefore calculated 
numerically, see below. The sum J2t=± refers to the two types of "almost critical" 
orbits with their crucial reflection point to the left (e = +) or to the right (e = — ) 
of the critical point A respectively. The Morse index o t is again defined by o t = 
// — |(sgn(/€ e ) — I) and U e is an abbreviation for U e = 7^7- Note that it is in general 
not allowed to replace the error function in eq. (10) by its asymptotic approximation 
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since U e = is possible for certain orbits, e.g., for those orbits which start and end 
in a self-conjugate point, see also [9]. 

In order to test this formula numerically, the contribution d co (E) of non-periodic 
critical orbits to the level density has been calculated using eq. (3). Integrating 
this density gives the contribution N co (E) of critical orbits to the spectral staircase 
function N(E) = ^{n\E n < E}. The integration constant is chosen such that 
lim^oo N co (E) = 0. 

The energy spectrum of the truncated hyperbola billiard has been calculated 
between E = and E = 10 000 (in natural units % = 2m = 1) using the boundary 
element method. In this interval all 1851 eigenvalues have been computed with high 
accuracy. No degeneracies have been found. 

In fig. 3 the fluctuating part N(E) — \ — N(E) is evaluated at the eigenvalues 
E = E n} with 

47T 47T 

being the generalized Weyl's law for this billiard system. Here A and L denote the 
area and perimeter respectively of the truncated hyperbola billiard and C 
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see 



[9]. The unusual long-range modulation of the spectral fluctuation visible in this fig- 
ure can be explained by the contribution N co (E) of critical orbits. The contribution 
of all 5 families of orbits with less than 5 reflections altogether has been calculated 
numerically. The total contribution is also presented in fig. 3, showing a very good 
agreement over the whole examined energy range. Subtracting the "correction term" 
N co (E) } see fig. 4, yields spectral fluctuations typical for billiard systems and should 
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semiclassically be explainable by Gutzwiller's periodic-orbit theory. 

The numerical results presented in this Letter clearly show that the contribution 
of non-periodic orbits to the integrated level density cannot be neglected for the 
truncated hyperbola billiard. The main reason for their numerical significance in 
the energy interval < E < fOOOO is their geometrical length of approximately - 
compared to the shortest periodic orbit with a length of 2. It is therefore expected 
that their effect becomes even bigger for larger values of a. We also predict an energy 
dependent critical orbit contribution to other quantities like the number variance, 
spectral rigidity, and the cosine-modulated heat kernel. 
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